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Abstract — The CALICE collaboration investigates different 
technology options for highly granular calorimeters for detec- 
tors at a future electron-positron collider. One of the devices 
constructed and tested by the collaboration is a lm 3 prototype 
for an imaging scintillator-steel sampling calorimeter for hadrons 
with analogue readout (AHCAL). The light from 7608 small 
scintillator cells is detected with silicon photomultipliers. The 
AHCAL has been successfully operated during electron and 
hadron test-beam measurements at DESY, CERN, and Fermilab 
since 2005. The collected data allow for evaluating the novel 
technologies employed. In addition, these data provide a valuable 
basis for validating pion cascade simulations. This paper presents 
the current status of comparisons between the AHCAL data 
and predictions from different Monte Carlo models implemented 
in Geant4. The comparisons cover the total visible energy, 
longitudinal and radial shower profiles, and the shower substruc- 
ture. Furthermore, this paper discusses a software compensation 
algorithm for improving the energy resolution of the AHCAL for 
single pions. 

Index Terms — GEANT4 validation, hadron calorimetry, silicon 
photomultipliers, CALICE, International Linear Collider. 



I. Introduction 

THE CALICE collaboration investigates new imaging 
calorimeter technologies for precision measurements at 
a future electron-positron linear collider. These calorimeters 
feature a fine granularity in both longitudinal and transverse 
direction. The high spacial resolution is needed to fulfil the 
shower separation requirement of Particle Flow reconstruction 
algorithms. The Particle Flow approach aims for a jet energy 
resolution of 3-4% at the proposed International Linear Col- 
lider ffl, 0, 0. 

CALICE has designed and constructed several prototypes 
for electromagnetic and hadron calorimeters. One of these 
devices is the prototype for an analogue hadron calorimeter 
(AHCAL). This prototype is a 1 m 3 scintillator-steel sampling 
calorimeter with 38 sensitive layers. The steel plates are 2 cm 
thick. Each scintillator layer is pieced together from separate 
tiles. The tiles are 5 mm thick and their lateral size varies 
between 3 x 3 cm 2 , 6 x 6 cm 2 , and 12 x 12 cm 2 . Figure [l] 
illustrates the layout of the sensitive layers. A wavelength- 
shifting fibre collects the light from each tile and guides it to 
a silicon photomultiplier (SIPM, c.f. (4)). The SIPM is directly 
attached to the tile. With a total of 7608 readout channels, the 
AHCAL prototype represents the first large-scale application 
of SIPMs. The detector is described in more detail in l5l . 
Figure [2] presents a cascade from a lOGeV pion measured 
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in the AHCAL during test-beam operation. This measurement 
illustrates the imaging capabilities of the device. 

The response of the calorimeter cells is equalised and 
calibrated using broad beams of muons acting as minimum 
ionising particles, i.e. the energy measured in each cell is 
expressed in multiples of the energy deposited by a muon in 
this cell (MIP). To suppress noise, only cells with a deposited 
energy of at least 0.5 MIP are considered for analysis. After 
correcting for the SIPM non-linearity, the uncertainty on the 
calibrated signal is roughly 3%. The electromagnetic scale, i.e. 
the conversion factor from MIP to GeV for an electromagnetic 
cascade, is determined using electron or positron data from 
test-beam measurements. Data sets covering different en- 
ergy ranges (energy and momentum are used interchangeably 
throughout this paper) and acquired under different operation 
conditions (CERN 2007: 10-50 GeV positrons; Fermilab 2009: 
1-20 GeV electrons) yield conversion factors that agree within 
their uncertainties. In addition, these data confirm the linearity 
of the electron (and positron) response up to 30 GeV. Figure [3] 
shows the AHCAL response to electrons and positrons of 
different beam momenta. The observed non-linearity above 
30 GeV (up to 3% at 50 GeV) is of the same order as the 
calibration uncertainty and is attributed to remaining effects 
from the SIPM saturation. More on the performance of the 
calorimeter can be found in (6) and (7). 

A detailed model of the AHCAL has been implemented 
in Mokka (version 7.02). Mokka is a Geant4 based 
Monte Carlo application and is capable of simulating test- 
beam setups and full detector geometries. The agreement 
between the MIP/GeV factor and other observables extracted 
from electron simulations and test-beam data validates the 
simulation of the various detector characteristics 0. 

This paper presents comparisons between data (collected 
with the AHCAL at CERN in 2007) and GEANT4 simula- 
tions for negative pions between 8 GeV and 80 GeV beam 
momentum. Several GEANT4 physics lists for the simula- 
tion of hadron interactions with matter exist. A physics list 
combines different models that are valid for different en- 
ergy ranges. At the transition between models, one model 
is randomly chosen for each incoming particle. This paper 
concentrates on four physics lists: QGSP_BERT, FTFP_BERT, 
FTF_BIC (GEANT4, version 9.3) and CHIPS (Geant4, ver- 
sion 9.3.p01). Validation studies performed by LHC experi- 
ments favour the QGSP_BERT physics list. Both FTFP_BERT 
and FTF_BIC have been re-tuned recently and are promising 
alternatives for QGSP_BERT. The CHIPS physics list is still 
at an experimental stage. However, this list is particularly 
interesting because it uses only one model for all hadron 




Fig. 1. The 38 scintillator layers of the AHCAL prototype. Each layer 
is 5 mm thick and is pieced together from 216 tiles measuring 3 x 3 cm 2 , 
6 x 6cm 2 , or 12 x 12cm 2 . 




Fig. 2. Measurement of the energy deposited by a lOGeV negative pion in 
the AHCAL. The primary ionisation track, the location of the first inelastic 
pion-nucleus interaction and the development of the hadronic cascade are 
visible. 



energies. Table [I] summarises the models used by these physics 
lists for simulating pions of different energies. A more detailed 
description of the physics lists and the individual models is 
given in ifTOl . 

II. Calorimeter Response to Pions 

Figure [4] presents the ratio between the mean visible energy 
for negative pions from simulation and measurements for beam 
energies between 8GeV and 80GeV IfTOl . The predictions 
made by the CHIPS physics list show an energy independent 
overestimation of roughly 8%. This overestimation is expected, 
because the low energy neutron cross-sections are not yet 
properly implemented in CHIPS. For the composite physics 
lists, the agreement between predictions and data depends on 
the particle energy. Both the QGSP_BERT and FTFP_BERT 
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Fig. 3. Mean reconstructed electron (Fermilab) and positron (CERN) 
energies at different beam momenta. The data confirm a linear response 
(indicated by the dashed line) up to 30 GeV. The non-linearity above 30 GeV 
is of the same order as the calibration uncertainty (3% at 50 GeV). 

TABLE I 

Composition of Geant4 physics lists. 



Physics list 


Model (for tt ± ) 


Energy range 


QGSP_BERT 


Bertini cascade (BERT) 
Low-energy parametrisation (LEP) 
Quark-gluon string model (QGSP) 


< 9.9 GeV 

9.5 GeV - 25 GeV 

> 12 GeV 


FTFP_BERT 


Bertini cascade (BERT) 
Fritiof string model (FTFP) 


< 5 GeV 
> 4 GeV 


FTF_BIC 


Binary cascade (BIC) 
Fritiof string model (FTF) 


< 5 GeV 
> 4 GeV 


CHIPS 


Chiral-invariant phase space model 


> OGeV 



physics lists agree with data within 2% at 8 GeV and overesti- 
mate the detector response by 6% at 80 GeV. Between 8 GeV 
and 25 GeV, QGSP_BERT shows a dip, which is attributed 
to the use of the LEP model as stopgap in this energy range. 
The FTF_BIC physics list underestimates the visible energy by 
4% at 8 GeV (Binary cascade and FTF) and gives 4% too high 
energy deposition at 80 GeV (FTF dominates the description 
of the first hard interaction). 

III. Topology of Pion Cascades 

Exploiting the high granularity of the AHCAL allows for 
determining the position of the first hard interaction of hadrons 
in the calorimeter (TT). According to Monte Carlo studies, the 
uncertainty of the algorithm applied for this determination is 
±1 calorimeter layer (1 layer w 3 cm) in about 74% of the 
events. The information on the position of the first interaction 
is used to investigate the average longitudinal shower profile 
relative to this point. Figure [5] presents a longitudinal shower 
profile relative to the calorimeter front face (filled histogram) 
and a longitudinal shower profile realtive to the first interac- 
tion point (open histogram). The latter profile looks shorter 
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Fig. 4. Ratio between the mean visible energy in the AHCAL from 
simulations and from data for negative pions at different beam energies. For 
all physics lists except for CHIPS, this ratio depends on the pion energy. 
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Fig. 5. Longitudinal shower profile for negative pions (45 GeV) relative 
to the calorimeter front (filled histogram) and relative to the first inelastic 
interaction (open histogram). 



and much more similar to an electromagnetic shower profile 
because the fluctuations of the shower starting position are 
removed. 



Figure 6a shows the longitudinal shower profile relative to 
the position of the first hard interaction for negative 1 8 GeV 
pions in the AHCAL ifTOl . Data are displayed as black points 
on top of the filled histogram from simulation (FTF_BIC). 
The general shapes of these profiles agree. In addition to the 
profiles, Fig.[6a]presents the breakdown of the energy contribu- 
tions from various particles to the shower (electrons, positrons, 
protons, and mesons). This additional information helps to 
discuss which physics process contributes most in which phase 
of the shower development. The overall longitudinal shower 
shape is dominated by electromagnetic processes (electron and 
positron contribution to the entire profile). The hadrons only 
contribute significantly in the very first layers after the first 



hard interaction. Figure [6b] shows the ratio between the mean 
of the longitudinal shower profiles (mean shower depth, centre 
of gravity along beam axis) from simulation and data. CHIPS 
overestimates the shower depth at energies below 80 GeV. 



The other physics lists underestimate the shower depth for 
all energies presented. 



Figure [7a] presents the radial shower profile for negative 
pions (18 GeV) in the AHCAL for data (black points) and 
simulation (FTF_BIC physics list, filled histogram). The over- 
all shapes of the profiles agree. Figure [7b] shows the ratio 
between the mean of the radial shower profiles (mean shower 
radius) from simulation and data. Except for QGSP_BERT 
below 10 GeV, all physics lists predict narrower showers than 
observed in data. CHIPS gives the best prediction of the 
shower radius. 

IV. Track Multiplicity 

The high granularity of the CALICE AHCAL provides 
the capability for identifying track segments from secondary 
hadrons produced within hadron showers lfT2l . A simple track- 
ing algorithm allows for finding tracks created by minimum 
ionising particles in the cascade. The algorithm relies on iso- 
lated hits and works on a layer-by-layer basis. The maximum 
angle 6 between the beam direction and the reconstructed 
tracks is intrinsically limited by the algorithm. For tracks in 
the 3 x 3 cm 2 tiles, the limit is #3 X 3 < 58°. For the larger tiles, 
the limit is correspondingly #6x6 < 72° and #12x12 < 81°. 
The fake track rate of this algorithm is on the level of few per 
mille. 

Figure [8] shows a typical hadron cascade from a 20 GeV 
negative pion in the AHCAL. Minimum ionising track seg- 
ments of the incoming pion and of secondary particles are 
identified and highlighted in red. The track multiplicity is 
influenced by the shower topology and especially by the 
number of secondaries created. The average track multiplicity 
is shown as a function of the beam energy for data and various 
simulations in Fig. [9] All physics lists predict too small track 
multiplicities. The physics list predicting values closest to data 
is QGSP_BERT. 

V. Performance and Software Compensation 

Hadron cascades consist of an electromagnetic and a 
hadronic component. The electromagnetic part originates from 
neutral pions and eta mesons created in inelastic hadron- 
nucleus collisions. These mesons decay practically instanta- 
neously into two photons, which give rise to electromagnetic 
cascades. Other hadrons involved in the cascade deposit part 
of their energy as invisible energy in the form of nuclear 
binding energy, nuclear recoil, or neutrinos. The AHCAL is a 
non-compensating calorimeter, i.e. the AHCAL response to 
the hadronic component of a pion cascade is smaller than 
the response to the electromagnetic component. The average 
electromagnetic fraction of a hadronic cascade increases with 
the momentum of the incident particle. The electromagnetic 
fraction and the invisible energy fluctuate strongly from one 
event to another. These two effects lead to a non-linear detector 
response to pions and degrade the detector resolution. Software 
compensation algorithms can be applied to correct for these 
effects. 

Hadron showers with higher energy density (shower energy 
divided by shower volume) tend to have a larger reconstructed 
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Fig. 6. (a) Longitudinal shower profile relative to the first interaction for negative pions (18 GeV) in the AHCAL for data (black points) and simulation 
(FTF_BIC physics list, filled histogram). The error bars include only the statistical uncertainty and the uncertainty introduced by the determination of the first 
inelastic scattering. The breakdown of the energy contribution from various particles in the shower (electrons, positrons, protons, and mesons) is shown, (b) 
Ratio between the mean of the longitudinal shower profiles (centre of gravity along beam axis) from simulation and data. CHIPS overestimates the shower 
depth, while the other physics lists behave opposite. 
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shower profile for negative pions (18 GeV) in the AHCAL for data (black points) and simulation (FTF_BIC physics list, filled histogram), 
the mean of the radial shower profiles (mean shower radius) from simulation and data. CHIPS gives the best prediction of the shower 



energy than showers with a smaller energy density for the 
same particle momentum. A clustering algorithm allows to 
select only energy depositions belonging to the actual hadron 
cascade and to exclude noise signals. Because of the high 
energy density in purely electromagnetic showers, a high 
energy density in a hadron shower is attributed to a large 
electromagnetic fraction, i.e. the cluster energy density gives 
an estimate of the electromagnetic fraction inside a cascade. 

Weights for the visible energy deposited by hadrons are 
extracted from Monte Carlo simulations based on the cluster 
density and the cluster energy. Applying individual energy 
weights (depending on both the intially measured cluster 
energy and the cluster density) to each pion event in the 
detector yields a significant improvement of the pion lin- 



normalised to the response to 15 GeV pions. The weights are 
extracted from simulations using the FTF_BIC physics list. 



earity and energy resolution. Figure 10 shows the linearity 
before (filled circles) and after (open circles) applying this 
software compensation technique. The AHCAL response is 



Figure 11 shows the AHCAL energy resolution for single 
pions before (filled circles) and after (open circles) applying 
software compensation. The compensation algorithm reduces 
the stochastic term (a/y/E) of the resolution by about 15%. 
Further information on the software compensation technique 
presented in this paper is given in fT3ll . Studies of more 
sophisitcated software compensation algorithms are ongoing. 

VI. Conclusions 

The CALICE collaboration has successfully built and op- 
erated an imaging calorimeter with analogue readout for 
hadrons. The pion data recorded with this detector allow for 
in-depth studies of the properties of hadronic cascades and the 
validation of Monte Carlo models with an unprecedented level 
of detail. With this highly granular calorimeter, it is possible 



to to determine the position of the first hard interaction, to 
investigate detailed shower profiles, and to measure the track 
multiplicity. The negative pion data presented in this paper 
cover the energy range from 8 GeV to 80 GeV. The analysis 
of pion data for lower energies down to 1 GeV is ongoing. 

The QGSP_BERT physics list shows a strong effect from 
the transition to the LEP parameterisation in the energy region 
between 9 GeV and 25 GeV. The effects from model transitions 
are reduced in the FTF-based physics lists (FTFP_BERT 
and FTF_BIC). For the observables presented in this paper, 
these physics lists yield the best overall performance with 
predictions agreeing with data within 5-10%. The CHIPS 
physics list is promising since it has no transition between 
different models, but this model is still in a testing phase and 
requires further development. 

The application of a simple software compensation algo- 
rithm to test-beam data improves the AHCAL energy resolu- 
tion for single pions by 15%. 
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Fig. 8. Energy deposited by a 20 GeV pion shower in the AHCAL. The 
signals identified as track segments are highlighted in red. 
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Fig. 9. Average track multiplicity at different beam energies extracted from 
data and predicted by several physics lists. 
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Fig. 10. AHCAL linearity extracted from test-beam data with a single energy 
weight (filled circles) and energy weights based on the cluster energy density 
(open circles). 
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Fig. 11. AHCAL energy resolution extracted from test-beam data with a 
single energy weight (filled circles) and energy weights based on the cluster 
energy density (open circles). 



